Spin Current Generation and Detection in the Presence of AC Gate 



o 
o 

(N 

> 
O 



A. G. Mal'shukov 1 , C. S. Tang 2 , C. S. Chu 3 and K. A. Chao 4 
1 Institute of Spectroscopy, Russian Academy of Science, 142190, Troitsk, Moscow oblast, Russia 
2 Physics Division, National Center for Theoretical Sciences, P.O. Box 2-131, Hsinchu 30013, Taiwan 
3 Department of Electrophysics, National Chiao-Tung University, Hsinchu 30010, Taiwan 
4 Solid State Theory Division, Department of Physics, Lund University, S-22362 Lund, Sweden 

We predict that in a narrow gap III-V semiconductor quantum well or quantum wire, an observable 
electron spin current can be generated with a time dependent gate to modify the Rashba spin-orbit 
coupling constant. Methods to rectify the so generated AC current are discussed. An all-electric 
method of spin current detection is suggested, which measures the voltage on the gate in the vicinity 
of a 2D electron gas carrying a time dependent spin current. Both the generation and detection do 
not involve any optical or magnetic mediator. 

PACS numbers: 71.70.Ej, 72.25.Dc, 73.63.-b 
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One key issue in spintronics based on semiconductor is 
the efficient control of the spin degrees of freedom. Datta 
and Das 1 suggested the use of gate voltage to control the 
strength of Rashba spin-orbit interaction (SOI) 2 which 
is strong in narrow gap semiconductor heterostructures. 
In InAs-based quantum wells a variation of 50% of the 
SOI coupling constant was observed experimentally. 3 4 
Consequently, much interest has been attracted to the 
realization of spin polarized transistors and other devices 
based on using electric gate to control the spin dependent 
transport £ 

Besides using a static gate to control the SOI strength 
and so control the stationary spin transport, new phys- 
ical phenomena can be observed in time dependent spin 
transport under the influence of a fast varying gate volt- 
age. Along this line, in this article we will consider a 
mechanism of AC spin current generation using time de- 
pendent gate. This mechanism employs a simple fact that 
the time variation of Rashba SOI creates a force which 
acts on opposite spin electrons in opposite directions. In- 
versely, when a gate is coupled to a nearby electron gas, 
the spin current in this electron gas also induces a vari- 
ation of the gate voltage, and hence affects the electric 
current in the gate circuit. We will use a simple model to 
clarify the principle of such a new detection mechanism 
without any optical or magnetic mediator. The systems 
to be studied will be ID electron gas in a semiconductor 
quantum wire (QWR) and 2D electron gas in a semicon- 
ductor quantum well (QW). 

We consider a model in which the Rashba SOI is de- 
scribed by the time dependent Hamiltonian H so (t) — 
ha(t)(k xP)-s, where k is the wave vector of an electron, 
Hs is the spin operator, and v is the unit vector. For a 
QWR v is perpendicular to the wire axis, and for a QW 
perpendicular to the interfaces. The time dependence of 
the coupling parameter a(t) is caused by a time depen- 
dent gate£ To explain clearly the physical mechanisms 
leading to the spin current generation, we will first con- 
sider the ID electron gas in a QWR, and assume a(t) to 
be a constant a for t<0, and a(i)=0 for t>0. For the ID 
system we choose the x direction as the QWR axis and y 
axis parallel to v, to write the SOI coupling in the form 



H ao (t)=ha(t)k x s z . For t<0 the spin degeneracy of con- 
duction electrons is lifted by SOI, producing a splitting 
A=Hak x between s z —l/2 and s z —— 1/2 bands, as shown 
in Fig. 1 by solid curves together with the Fermi energy 
Ep. The spin current in this state is zero, as it should 
be under thermal equilibrium. 

Indeed, the spin current is defined as I s (t)=I^(t) — 
Il(t), where 2j(i) [or is the partial current asso- 

ciated with the spin projections s z =l/2 (or s z =— 1/2). 
Hence, 
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where L is the length of the QWR. Taking the momentum 
derivative of the Hamiltonian, we obtain the velocity as 



The 



v vi {k x ) = hk x /m* ±a{t)/2. (2) 
spin current is then readily obtained as 

I s (t) = {hn/Am*){Kk^ - Kk{) + Ha(t)n/4, (3) 



where n is the ID electron density, and fej (or fejj is the 
average momentum in the f-spin (or |-spin) band. 

For a parabolic band fifc-f =— m*cv/2 and hki—m*a/2. 
Although Hkij — hki gives a finite contribution to / s (t) in 
P)l. for t<0 where a(t)=a, this contribution is compen- 
sated by the contribution Han/ 4 due to the SOI. Hence, 
the total spin current l s (t)=0 for i<0. However, when 
the SOI is switched off at i=0, a(t)=0 and so the spin 
current is finite, because the average electron momenta 
retain the same as they were at t<0. As time goes on, the 
electron momenta relax with a relaxation time r. There- 
fore, J a (i)=— (fam/4) exp(-i/r) for t > 0. 

It is instructive to make a Fourier transform of I s {t) 
to obtain a Drude-like expression 
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Since the units of our spin current is h/2, the above ex- 
pression is a complete analogy to the electric conductiv- 
ity. Instead of an electric driving force eE, here we have 
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an equivalent driving force (m* /2)[da(t) / dt], the Fourier 
component of which is (m*/2) zf2 a(fi). Under this driv- 
ing force we have the classic equation of motion 



dv^i _^m* da(t) 



dt 



(5) 



This force acts in opposite directions on electrons with 
opposite spin projections. When such a force creates a 
spin current, it does not induce an electric current. 

The above conclusion of spin current generation can 
be demonstrated with a rigorous linear response analysis, 
which will be performed on a 2D electron gas (2DEG). 
The simple Drude expression Q will then appear as a 
general result. Let the 2DEG be in the xy plane with the 
unit vector v along the z axis, which is the spin quanti- 
zation axis. We will use the equation of motion for the 
spin density operator to generalize the ID expressions 
|2J) for the spin current. For a homogeneous system the 
spin current density operators can be expressed in terms 
of the electron creation operator ct and destruction op- 
erator eg , where 7 labels the spin projection onto the 
z-axis. This current is then derived as 
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where the superscript i—x,y,z specifies the direction of 
spin polarization, and the subscript j=x,y refers to the 
direction of the spin current flow. The partial current 



When an AC bias with frequency is applied to the front 
or the back gate of a 2DEG, 3,4 the Rashba coupling con- 
stant contains two terms a(t)=cto+5a(t) 1 where ao is 
constant in time and 5a(t)=6a e . We assume that the 
only effect of the AC bias is to add a time dependent com- 
ponent to the SOI coupling constant, although in practice 
it is not simple to avoid the bias effect on the electron 
density** The SOI Hamiltonian is separated correspond- 
ingly into two parts H S0 (t)=Hg +H' S0 (t). The time inde- 
pendent part H® does not produce a net spin current in 
the thermodynamically equilibrium state. However, as 
pointed out in the above analysis on the 1DEG system, 
the time dependent H so (t) can give rise to a spin current. 

We will incorporate H° into our unperturbed Hamil- 
tonian and treat H' so (t) within the linear response regime. 
The so-generated AC spin current (jHt)) has the form 



(j;(t)) = \ \ dt'{[H' so {t'),j]{t)}) 



+ e llz HSa(t) n/4. 



(10) 



In the above equation the first term can be written in the 
form 6a(t)lZ l j(n). For zero temperature and with £!>0, 
the response function can be represented as the 

Fourier transform of the correlator 
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is the ordinary kinematic term, and 
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is the contribution of SOIj£ Here e uz denotes the Levy- 
Civita symbol. The SOI induced current resembles the 
diamagnetic current of electrons under the action of an 
external electromagnetic vector potential. 

We note that the SOI Hamiltonian can be conveniently 
written in terms of the kinematic current as 



H so (t) = [m*a(t)/h] (Jy - jy) . 
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FIG. 1: The dashed curve is the electron energy band without 
SOI. The SOI splits the energy band into the j-spin and the 
j-spin band, as shown by the solid curves, with corresponding 



average wave vector fcj and k± 
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where h^=kxv. In the above equation, the bar over the 
product of two one-particle Green functions means an 
ensemble average over impurity positions. 

We will use the standard perturbation theor^ to cal- 
culate this ensemble average, which is valid when the 
elastic scattering time r due to impurities is sufficiently 
long such that EpT^>h. We will assume that the elec- 
tron Fermi energy Ep is much larger than both Ml and 
haoh^. To the first order approximation, we neglect the 
weak localization corrections to the correlator , since 
these corrections simply renormalize the spin diffusion 
constant. 8 Consequently, the configuration average of the 
pair product of Green functions is expressed in the so- 
called ladder seriesi* We found that since he = 
many of such ladder diagrams vanish after angular inte- 
gration in l|llf> . similar to suppression of ladders in the 
electric current driven by the vector-potential^ At the 
same time, some of nondiagonal on spin indices diagrams 
do not turn to after the angular integration. Employ- 
ing the analysis of similar diagrams done in it can be 
shown that they cancel each others Hence, the configu- 
ration average in decouples into a product of average 
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Green functions and Eq. becomes 

^ (n) = - i ^E £ ' nz E fc ^« ( 12 ) 

l,n £ 
x /" < ^-Tv\ [ s l G{k,uj)s l G{k,Lo + n) , 

where G(fc, w) is the average Green's function which con- 
tains fully the effect of H® - This function is represented 
by the 2x2 matrix 



G{k,ui) = uj — E^/h — aoh^ ■ s + i rsgn(w) 



(13) 



where r=l/2r, and E^ i s defined with respect to Ep. 
Substituting (|13(l into Ijl2|l . and then into l|l()|l. we obtain 
the spin current 
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It is important to point out that the spin density under 
the gate area is zero. This is the reason why even in 
a 2DEG the D'yakonov-Perel spin relaxationiS does not 
appear in (|14fl for the generated spin current, although 
this spin current is determined by the response function 
(|12fl which involves spin degrees of freedom. Hence, in 
the homogeneous system with zero spin density, only elec- 
tron momentum relaxation occurs in the process of spin 
current generation by a time dependent gate. 

Unlike the spin current in a ID system, in a 2DEG 
the current given by (|14(l has no specific direction. To 
clarify the spatial distribution of the spin flux induced 
by an AC gate, let us take the chiral component J c hii(t) 
of the spin current 



J C KAt) = [(J v x (t))-(jy(t))}/2. 



(15) 



It is easily to see that this chiral projection has the same 
form as the expression for a ID system, if n represents 
the electron density of the 2DEG. In Fig. 2 we illustrate 
the spin current distribution for a circular gate which is 
marked as the gray area. The spin polarization at any 
point under the gate has two components parallel to the 
2DEG. For any direction specified by the unit vector N, 
the two spin polarized flux with polarization directions 
parallel and antiparallel to N will oscillate out of phase by 
the amount of tt along the direction perpendicular to N. 
Such out of phase oscillation is schematically plotted in 
Fig. 2. The amplitude of the spin density flow in each of 
the opposite directions, as marked by the dashed-line ar- 
rows, is just J'chir(i)- In the 2DEG outside the gate area, 
the spin current can be supported only by spin diffusion. 
Therefore the chiral AC spin polarization is accumulated 
in the vicinity of the circumference of the gate, and from 
where diffuses away from the gate area. It can also dif- 
fuse under the gate. For small gates such back diffusion 
can diminish the efficiency of the spin generation. On 
the other hand, for large gates with the size larger than 




FIG. 2: Distribution of spin currents induced by a time de- 
pendent circular gate which is marked as the gray region. 
Under the gate, electrons with opposite spins (solid arrows) 
move in opposite directions indicated by the dashed-line ar- 
rows. Arrows outside the gate area show the accumulated 
spin polarization during a half period of AC gate voltage os- 
cillation. 



the spin diffusion length the diffusion counterflow does 
not reduce much the total spin current. 

The so-generated current amplitude can be easily 
estimated. With <Sa=3xl0 6 cm/s, 4 for fi=27rxl0 9 
s _1 , n=10 12 cm -2 , and r=l ps, from l|14fl we derive 
(2e/ft)(J J i (fi))~10^ 3 Amp/cm. This AC spin current 
can be detected by various methods. For example, if 
holes can tunnel into the neighborhood of the gate edge, 
their recombination with spin polarized electrons will 
produce the emission of circular polarized light". 

However, here we will discuss a new method of direct 
electric detection of the DC or the AC spin current. This 
method is based on a simple fact that the Rashba SOI 
couples the spin current to the gate voltage. We have 
shown in our above analysis that due to this coupling, 
spin current can be induced by a time dependent gate 
voltage. In this case the voltage variation plays the role 
of a source which drives electrons out of thermodynamic 
equilibrium, and the spin current is the linear response to 
this perturbation. The reverse process is to create a spin 
current in a 2DEG by some source, and so inducing a 
voltage shift in a nearby gate. This is also possible to re- 
alize. We thus consider a model where the SOI constant 
a(U) is a function of the gate voltage U(t)—Uo+V(t). 
Uq is the static equilibrium value in the absence of a spin 
current, while V(t) is a dynamic variable. The mean 
value (V) of V{t) has to be calculated as a linear re- 
sponse to the perturbation associated with the presence 
of the spin polarization flow. The explicit form of this 
perturbation can be obtained by averaging the Hamilto- 
nian of the system over an electronic state with the given 
time dependent spin current. 

Let (■■■).] be such type of average. To the lowest order 
with respect to SOI, the coupling of the gate voltage to 
the spin current is thus determined by the average of the 
Rashba interaction in 10 with a=a(U). The coupling 
between the gate voltage U (t) and the spin current is 
via the kinetic current Jj. To derive the coupling Hamil- 
tonian H int , we use © to express Jj in terms of J]j, 
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FIG. 3: Schematic illustration of spin current detection. AC 
spin current flows from the right to the left under the gate 
with spin polarized as shown by arrows. V denotes the volt- 
meter and Z is the outer circuit impedance. 



and expand a(U)=a(Uo)+a'V(t) for small V(t). 
coupling Hamiltonian is then derived from J^J) as 



Hi, 
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The charging of the gate Q=CV is related to the gate 
capacitance. Hence, (TTfi)) can be expressed in the conve- 
nient form H int =Q£, where 



{J?). 



(17) 



is the effective electromotive force. 

To illustrate our proposed method of direct electric de- 
tection, let us consider a circuit connected to the gate. 
The principal scheme of the spin current detection is 
shown in Fig. [3] In it, an additional back gate can be 
utilized to tunc the electron density (not shown) . The cir- 
cuit is characterized by a frequency dependent impedance 
Z(£l). The voltage induced on the gate by the electro- 
motive force Ijl7|l is then easily obtained as 



(V)=£ 



i + incz(n) 



(18) 



When the spin current frequency is in resonance with the 
circuit eigenmode, the gate voltage becomes very large. 
In the limit of high impedance (open circuit), (V)=£. 
Using the spin current (2e/fi.)(j7J(f2)}~10~ 3 Amp/cm de- 
rived above, and the fact that (Jl)j — A (J^(O)), where 
A is the area under the gate, let us estimate the electro- 
motive force induced in a probe gate by this spin current 
generated by a nearby source gate. For the reasonable 
parameter values cv'=3xl0 7 cm/Vs^ m*=0.03 m e , and 
C~K6 A/l with k=10 and Z=1CT 5 cm, from (|14"}| and 
JTHJ we obtain £~1CT 5 Volts. 



The generated AC spin current can be rectified with 
various methods. For example, one can use a shutter gate 
which is 7r/2 phase shifted with respect to the generation 
gate. The shutter gate can be placed in the neighbor- 
hood of the generation gate or between two such gates. 
The evaluation of the rectifying efficiency of such a setup 
requires a thorough analysis of spin relaxation and dif- 
fusion processes caused by the spin accumulation during 
the shutter cycle. 



We would like to add one relevant information which 
we became aware of after we have completed this paper. 
The preprint of Governale et al. on the quantum spin 
pumping in a ID wire is also based on the idea of cre- 
ating spin current via time dependent gate^ However, 
our results involving dissipative transport in 2DEG and 
1DEG can not be compared directly with those in Ref.lT^. 
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